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An efficient multi-mode quantum memory is a crucial resource for long-distance quantum com- 
munication based on quantum repeaters. We propose a quantum memory based on spectral shaping 
of an inhomogeneously broadened optical transition into an atomic frequency comb (AFC). The 
spectral width of the AFC allows efficient storage of multiple temporal modes, without the need 
to increase the absorption depth of the storage material, in contrast to previously known quan- 
tum memories. Efficient readout is possible thanks to rephasing of the atomic dipoles due to the 
AFC structure. Long-time storage and on-demand readout is achieved by use of spin-states in a 
lambda-type configuration. We show that an AFC quantum memory realized in solids doped with 
rare-earth-metal ions could store hundreds of modes or more with close to unit efficiency, for material 
parameters achievable today. 



I. INTRODUCTION 

The distribution of entanglement between remote lo- 
cations is critical for future long-distance quantum net- 
works and extended tests of quantum non-locality. It 
is likely to rely on quantum repeaters P, [2], which re- 
quire quantum memories that can store entanglement 
between distant network nodes [1, Recent experi- 
mental achievements in quantum state storage d, H, 0, [1| 
demonstrate that currently investigated QMs can store a 
single mode. Yet, long-distance quantum repeaters hav- 
ing QMs only capable of storing one mode would only 
generate very limited entanglement generation rates [9|. 
To achieve useful rates some way of multiplexing the 
QM will be required [l^, [HI ■ By using time [ll| , spa- 
tial jl^j [13, I O'" frequency multiplexing to store sin- 
gle photons in many modes N, the entanglement gen- 
eration rate can be increased by a corresponding factor 
N [llj. Here we consider multi-mode QMs capable of 
storing N temporally distinguishable modes, which is 
a natural form of sending information also used in to- 
day's telecommunications networks. Time multiplexing 
is extremely challenging using current QM protocols such 
as stopped light based on Electromagnetically Induced 
Transparency (EIT) photon echoes based on Con- 
trolled Reversible Inhomogencous Broadening (CRIB) 
[TgI or off-resonant Raman interactions [l^. For a 
single mode, iV = 1, the efficiency as a function of op- 
tical absorption depth d is similar for EIT, CRIB, and 
Raman-type memories [H, [l^ . For multi-mode storage 
using EIT and Raman, however, N scales as \Q [20l.[2ll|. 
which severely limits the number of potential modes in a 
practical experiment, while CRIB offers a better scaling 
since iV ~ d [m, iH. Still, even using CRIB, efficient 
multi-mode storage requires extreme absorption depths 
that are currently out of reach. 

In this article we present a QM, where the number 
of potential modes is essentially independent of the op- 
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tical depth. As we will show, efficient storage of 100s 
of modes or more is then feasible with realistic material 
parameters. The QM is based on absorption of light by 
an ensemble of atoms on an optical transition which is 
inhomogeneously broadened. Ensembles are in general 
very attractive as QMs due to strongcoUective enhance- 
ment of the light- matter coupling Storage of sin- 
gle photons using EIT (stopped light) has been demon- 
strated with cold alkali atoms [6,, 3] which can be treated 
as homogeneous ensembles of identical atoms. Here we 
will instead consider ensembles of rare-earth-metal (RE) 
ions in solids, which are inhomogeneously broadened. 
Stopped fight with storage times up to 1 s has been 
demonstrated in RE-doped solids [22, [23| , where an ap- 
proximate homogeneous ensemble was created by spec- 
trally isolating a narrow absorption peak through opti- 
cal pumping. In contrast to this approach, the quantum 
memory proposed here uses the inhomogeneous broaden- 
ing as a resource in order to achieve better multi-mode 
performance. For that one needs to coherently control 
the dephasing which is caused by the inhomogeneous fre- 
quency distribution of the atoms. The key feature of our 
proposal is to achieve this control by a specific shaping 
of this distribution into an atomic frequency comb. 

The article is organized in the following way. In Section 
II we give an overview of the proposal using a simplified 
physical picture. In Section III we show results from 
an analytical analysis of the physics of atomic frequency 
combs, which is further detailed in the Appendix. In 
Section IV the multi-mode storage capacity is discussed, 
and in Sec. V we discuss implementation of the quantum 
memory in rare-earth-ion-dopcd solids. Conclusions are 
given in Sec. VI. 



II. ATOMIC FREQUENCY COMBS 

The storage material is an ensemble of atoms with an 
excited state |e) optically connected to two lower states 
\g) and |s), see Fig. [TJ The states \g) and |s) are typ- 
ically ground-state hyperfine or Zeeman states. We as- 
sume that the optical transition \g) — \e) has a narrow 
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homogeneous linewidth 7/1 and a large inhomogeneous 
broadening Tin {Tm/lh ^1)- This is the case for RE- 
doped soHds that we will consider for physical implemen- 
tation in Sec. V. The \g) — |e) transition is spectrally 
shaped such that the atomic density function consists 
of a series of narrow peaks spanning a large frequency 
range (atomic frequency comb, or AFC). This can be 
done by frequency-selectively transferring atoms from \g) 
to a metastable state \aux), for instance a third hyper- 
fine state, through optical pumping techniques (see Sec. 
V). Note that the maximum spectral width of the AFC 
is limited by the level spacings between the ground state 
levels. 

We then consider a single-photon input field, having a 
spectral distribution 7p larger than the AFC peak sepa- 
ration A, but narrower than the total width of the AFC 
F, which is in resonance with the |g) — |e) transition. If 
the atomic density integrated over the photon bandwidth 
is high enough the photon can be totally absorbed by the 
AFC, although the spectral density of atoms is concen- 
trated to narrow peaks. This can be understood in terms 
of the Heisenberg energy-time uncertainty relation. For 
the time scale of the absorption, which is of the order of 
the input pulse duration r = l/7p, the optical transition 
will have an uncertainty of the order of 7p ^ A. This 
causes a spectral averaging of the sharp AFC structure 
into a smooth distribution, allowing for uniform absorp- 
tion over the photons bandwidth (see also Sec. III). 

After absorption (at t — 0) the photon is stored as a 
single excitation de-localized over all the atoms in the 
ensemble that are in resonance with the photon. The 
state is described by a collective Dicke state [24 1 
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FIG. 1: (Color online) The principles of the proposed atomic 
frequency comb (AFC) quantum memory, (a) An inhomoge- 
neously broadened optical transition \g) — \e) is shaped into 
an AFC by frequency-selective optical pumping to the \aux) 
level. The peaks in the AFC have width (FWHM) 7 and are 
separated by A, where we define the comb finesse as F = A/7, 
(b) The input mode is completely absorbed and coherently ex- 
cites the AFC modes, which will dephase and then rephase 
after a time 2n/A, resulting in a photon-echo type coherent 
emission. A pair of control fields on |e) — \s) allow for long- 
time storage as a collective spin wave in |s), and on-demand 
read-out after a storage time Ts- 
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where Zj is the position of atom j, fc is the wave- number 
of the light field (for simplicity, we only consider a single 
spatial mode defined by the direction of propagation of 
the input field) , Sj the detuning of the atom with respect 
to the laser frequency and the amplitudes Cj depend on 
the frequency and on the spatial position of the particular 
atom j. 

The collective state can be understood as a coherent 
excitation of a large number of AFC modes by a single 
photon. These modes are initially (at t = 0) in phase 
with respect to the spatial mode fc. But the collective 
state will rapidly dephase into a non-collective state that 
does not lead to a strong collective emission, since each 
term acquires a phase exp{i6jt) depending on the de- 
tuning 5j of each excited atom. If we consider an AFC 
having very sharp peaks, then the detunings Sj are ap- 
proximately a discrete set such that 6j — nijA, where 
rrij are integers. It follows that the collective state is re- 
established after a time 27r/A, which leads to a coherent 
photon-echo [1^ |2^ |23| type re-emission in the forward 
direction. The efficiency of this process (see Sec. Ill for 



details) can reach 54% in the forward direction (limited 
only by re-absorption) . But if the re-emission is forced to 
propagate in the backward direction, by a proper phase 
matching operation (see below), the process can reach 
100% efficiency. 

The process described so far only implements a QM 
with a fixed storage time. In order to allow for on- 
demand read-out of the stored field (which is a necessary 
requirement for use in quantum repeaters) and long-term 
storage, the single collective excitation in |e) is trans- 
ferred to a ground state spin level \s). This can be done 
by applying an optical control field on |e) — |s), for in- 
stance a short tt pulse. The excitation is thereafter stored 
as a collective spin wave, which is also the basis of stor- 
age in FIT, CRIB, and off-resonant Raman type mem- 
ories based on ensembles of lambda atoms [ij, [l^ Iv^ . 
The spin wave allows for long-lived storage since spin 
coherence lifetimes are generally longer than the optical 
coherence lifetimes. Moreover^ the spin transition does 
not have a comb structure [28||, which means that the 
evolution of the phases in Eq. ([T]) will be locked during 
the storage in the spin wave. If we assume that the spin 
transition is completely homogeneous the storage time is 
only limited by the coherence lifetime of the spin wave. If 
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there is inhomogeneous broadening, spin echo techniques 
can be used to exploit the homogeneous coherence hfe- 
time of the system (see Sec. V). To retrieve the field a 
second counter-propagating control field is applied after 
some time Tg . These operations will result in a backward 
propagating output mode after a total time Tg + 27r/A. 
As we will show below, the efficiency of the storage for 
backward retrieval can theoretically approach unity for 
high optical depths. 



III. ANALYTICAL TREATMENT 

We will now discuss some important features of the 
AFC and give quantitative formulas relating the input 
field to the output field as a function of the AFC param- 
eters. This can be used to calculate, for instance, storage 
efficiency. To this end we have used an approximate, lin- 
earized Maxwell-Bloch (MB) model. A linearized model 
can be used since the population in the excited state is 
negligible. For a single photon input field at most one 
atom in the ensemble is excited. For multi-mode stor- 
age, the number of totally excited atoms should be much 
less than the total number of atoms, which is the case 
for the number modes we consider and the RE-solids we 
discuss in Sec. V. The linearity of the analytical model 
also implies the equivalence between classical and single- 
photon dynamics ^17i . JS . i29i] . The important features 
of the model is discussed below, while the details of the 
calculations are given in the appendix. 

The matter part was assumed to be an ensemble of 
two- level atoms, including levels \g) and |e) in Fig. [U 
and the light field was treated as a one-dimensional spa- 
tial mode. This model can fully describe the mapping 
of light onto the atoms, the atomic evolution including 
the inhomogeneous dephasing and rephasing of the AFC 
modes, and the re-emission of light. The model did not 
include population relaxation (Ti) or decoherence (2^2 ). 
This is a valid approximation since the peaks in the AFC 
are assumed to be inhomogeneously broadened. Then all 
relevant dephasing times are much faster than Ti and 
T2 (as in the example given in Sec. V). In practice, 
however, the width of the individual peaks will be lim- 
ited by the homogeneous linewidth 7/j = l/(7rT2). The 
two counter-propagating control fields were modelled by 
transferring the forward atomic mode, associated with 
the forward propagating optical mode, to a backward 
atomic mode [2^. This achieves the phase matching op- 
eration that results in a backward propagating optical 
mode. Clearly this approach assumes that the transfer 
efficiency of the control pulses is unity and that the spin 
transition \g) — \s) is decoherence free. We thus only ana- 
lytically treat the physics of the atomic frequency comb, 
which is the novel feature of this protocol as compared 
to previous quantum memory proposals [3, [H, [13] ■ The 
transfer efficiency and spin decoherence could, however, 
be straightforwardly included with appropriate loss fac- 
tors. The equations of motion describing the light-matter 



dynamics is given in the appendix (see Eqs (|A1[) - (|A4[) '). 
In addition to the analytical calculations, we also show 
results from numerical solutions of the MB model. 

The atomic spectral distribution n(6) is a frequency 
comb described by a series of Gaussian functions charac- 
terized by the peak separation A and peak width 7, with 
an overall comb width of F 



n(<5) cx e-*'/(2i^') ^ g-(^-,A)V(27^) (2) 

j = -oo 

The full- width at half-maximum (FWHM) peak width is 
then 7 — v81n27, and we define the finesse F of the 
comb as F = A/7. We are interested in the regime 
where we have a significant number of well-separated 
peaks within the overall characteristic width of the AFC. 
Moreover, we assume that the spectrum of the photon is 
narrower than the width of the AFC, but wider than the 
peak separation, i.e. F 7^ ^ A ^ 7. 

In order to have an efficient memory, the probability of 
absorption of the input mode must be sufficiently high. 
But to create a high finesse AFC atoms are removed, by 
optical pumping, from the natural inhomogeneous spec- 
trum. This preparation step lowers the total density 
of atoms and hence the probability of absorption. The 
analytical calculations show that the amplitude of the 
forward propagating mode Ef{z,t) after the sample of 
length L is given by (see Eqs ^M^-^ME^) 

Ef{L,t) = EfiO,t)e-'^ = EfiO,t)e-^, (3) 

where the effective absorption depth d — aL can be ex- 
pressed as 

d^i.r^^i, (4) 

F V 41n2 F ^ ' 

The comb absorption is thus reduced by a factor of 
1/F, which is to be expected since the total density of 
atoms in the comb (integrated over the frequency band- 
width of the photon) is reduced by a factor of 1/F in the 
preparation step. But if the peak absorption depth d of 
the comb is high enough, the total density of atoms in 
the comb can still be sufficiently high to absorb the input 
mode with close to 100% probability. 

An interesting feature of the AFC is that it can ab- 
sorb uniformly over the frequency spectrum of the in- 
put mode although the atomic frequency comb has large 
" holes" . This can be understood by considering that the 
absorption is an event well localized in time, of the order 
T = l/7p. In that case the light-atom system reacts with 
a Fourier- limited resolution 7p, which effectively smears 
out the structure of the AFC to an essentially fiat distri- 
bution (see discussion preceeding Eq. (jAlip V In other 
words the peak separation A must be smaller than the 
photon bandwidth 7p. This also implies that Eq. ([3]) is 
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only valid for times t <C Stt/A, i.e. before the appear- 
ance of the first emission due to rephasing of modes in 
the comb. 

By applying the two counter-propagating control 
fields, the spatial phase pattern of the atomic polariza- 
tion changes such that it couples with the backward prop- 
agating field mode [l^ [l^ |23| . Under the assumption of 
unit transfer efficiency and no phase decoherence during 
the transfer to or storage in |s), the atomic phases will 
evolve to a collective state coupled to the backward prop- 
agating mode. This leads to emission in the backward 
mode. In the appendix (sec Eq. (jAlSp ) we show that the 
output field amplitude is a coherent sum over the ampli- 
tudes corresponding to all possible paths of the photon 
in the medium, where each path is characterized by a 
particular turning point in the medium. Even though 
there is absorption also along the path of the backward 
mode, a constructive interference effect between all these 
possible paths leads to the output amplitude potentially 
being arbitrarily close in magnitude to the input ampli- 
tude for large enough d. This interference effect is also 
present in the CRIB quantum memory scheme for back- 
ward retrieval [29]. The analytical calculation leads to 
the following relation between the forward input field 
Ein — Ef{z = 0,t — 0) and output backward field 
Font ^Ebiz = 0,t^ 2tt/A) (see Eq. ^MB), 

Eout = -£;„e^^2'^^(l - e-^)e-^^. (5) 

The first factor represents a global phase shift that de- 
pends on the relative detuning Aq between the carrier 
frequency of the light and the centre of the AFC. The 
second factor describes the coupling between the input 
and output fields, which increases with increasing effec- 
tive absorption depth, Eq. (U), whereas the third factor 
takes into account dephasing due to the finite width of 
the peaks. 

The dephasing factor can be understood by considering 
the dephasing during the storage in the optically excited 
state |e). In the qualitative description in Sec. II we 
assumed that each peak was a delta function, which leads 
to complete rephasing of the collective state Eq. ([T]). 
Clearly finite Gaussian peaks will not lead to complete 
rephasing, which lowers the probability of re-emission. 
The decay of the re-emitted field amplitude is given by 
the Fourier-transform of a peak in the comb (see Eq. (jA6|) 
and discussion preceding Eq. (|A16[) ). in analogy with 
the role of the initial peak shape in the CRIB quantum 
memory protocol [1^ [s^l • In the case of Gaussian peaks, 
the decay is given by e~* ^ where t = 27r/A is the 
storage time in the excited state. This factor can be 
expressed in terms of the finesse F, which leads to the 
last term in Eq. ([5]) . To reduce this dephasing one should 
thus maximize the comb finesse F. 

The efficiency, defined as 77 = \Fout\^ /\Ein\'^ , is plotted 
as a function of optical depth for different comb finesses 
in Fig. [21 For a moderate finesse of = 4, an efficiency 
greater than 50% can be achieved for d — 10. For a high 
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FIG. 2: (Color online) The storage efficiency is shown as a 
function of peak absorption depth d for AFC finesse values 
of -F = 4, 6 and 10. The solid lines represent the efliciency 
calculated using the analytical result, Eq. ((5)l, whereas the 
symbols represent the numerical calculation. As an example, 
for d = 40 and F = 10, an efficiency of 90% can be achieved. 
Even for an absorption depth of 25 and a modest finesse of 6, 
the efficiency reaches 80%. 
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FIG. 3: (Color online) Efficiency as a function of comb finesse 
F for peak absorption depths d = 5 (black solid line), d = 10 
(red dashed line), d = 20 (green dotted line) and d = 40 (blue 
dashed-dotted fine). The curves were calculated using Eq. 



finesse, F > 10, the efficiency can be close to 100% for 
large enough d. From an experimental point of view d is 
often a fixed parameter. The storage efficiency can then 
be maximized for the given material by optimizing F us- 
ing Eq. (O (see Fig. [S]). This optimum in the choice 
of F is due to a trade-off between absorption probabil- 
ity (favoring small F) and low dephasing during storage 
(favoring large F). 
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We will also briefly discuss the possibility of re- 
emission in the forward direction. This will occur if no 
control fields are applied or if the control fields are ap- 
plied co-propagating. As is shown in the Appendix (sec 
Eq. (jA19|) ). this re-emission is limited in efficiency to 
54% for an optimal optical depth of d — 2. The optimum 
represents a trade-off between strong light-matter cou- 
pling (favoring high optical depth) and re-absorption of 
the re-emitted light as it propagates through the sample 
(favoring low optical depth). 



IV. MULTI-MODE STORAGE CAPACITY 

The excellent multi-mode properties of the AFC based 
QM can be understood qualitatively as follows. We imag- 
ine that the input is a series of temporally distinguish- 
able modes containing, with a certain probability, a single 
photon in each mode. The storage efficiency is uniform 
over all the modes, since each mode spends the same 
time in the memory. Note that the effect of spin deco- 
herence is also the same on each mode. If the control 
fields are applied as late as possible, the total duration of 
the input can be close to T = 27r/A, the time it takes to 
rephasc the first mode that was absorbed. The shortest 
duration r of one mode, however, is only limited by the 
total frequency bandwidth r ~ l/(A"pA), where Np is 
the number of peaks in the comb. The number of modes 
N — T/t ^ Np is then proportional to Np, such that 
adding more peaks in the AFC allows one to store more 
modes. Thus, in materials with inherently large inho- 
mogeneous broadenings (see Sec. V) one can increase 
the multi-mode capacity by making use of a wider range 
of the inhomogeneous broadening, without having to in- 
crease the peak absorption depth by, for instance, in- 
creasing the atomic density in the sample. 

We can make a more detailed calculation of the multi- 
mode capacity by assuming that the input modes have 
Gaussian temporal shapes, where each mode occupies a 
temporal slot of duration r. It is important that the 
modes do not overlap too much, which in the end would 
produce errors. The amount of overlap that can be ac- 
cepted depends on the particular application. In Ref. 
[Ill ] an analysis of the error due to overlap was carried out 
in the context of a quantum repeater protocol in which 
the AFC QM could be used. There it was shown that 
T should be chosen as r = 127r/r, where F is the total 
memory bandwidth (note that here F is expressed in ra- 
dians) , where the width F is taken to be that of a square 
distribution, as in Ref. This also guarantees that 

the spectral content of the input modes can be fitted in- 
side the AFC, which in turn results in negligible temporal 
distortion of the output modes. The width F can be re- 
lated to the number of AFC peaks A'^p = F/A. The total 
duration T of the input is T = 27r/A — Tq, where Tq 
is the time necessary for applying the two control fields. 
In the limit where Tq < 27r/A, we have T « 27r/A. 
The number of total temporal modes then amounts to 



N = T/tk Np/6. 

The number of peaks that can be created will depend 
strongly on the inhomogeneous Ti„ and homogeneous 
linewidths of the material. The width of each peak can- 
not be narrower than the homogeneous linewidth, i.e. 
7 > 7ft. Since an efficient memory (> 90%) requires a 
finesse F > 10, it follows that A > 10^ h- The num- 
ber of peaks will depend on the ratio of the maximum 
AFC bandwidth F to the peak separation A. The band- 
width will ultimately be limited by the inhomogeneous 
broadening of the material. In most materials, however, 
it will instead be limited by hyperfine or Zeeman transi- 
tion spacings in order to avoid exciting multiple transi- 
tions that would induce quantum beat effects. In Section 
V we will give a detailed example where the number of 
peaks is of the order of hundreds with known material 
parameters. 

The AFC QM dramatically reduces the necessary ab- 
sorption depth for multi-mode storage as compared to 
CRIB or FIT QMs, since the multi-mode capacity is in- 
dependent of the peak absorption depth of the material. 
For CRIB, d — 30N is required for storing N modes 
with an average efficiency of 90% [il|, so for instance 
N = 100 requires d = 3000. In the case of FIT even 
larger d would be necessary due to the worse scaling as 
mentioned in the introduction. The AFC only requires 
d = 40 to reach the same efficiency, independently of N. 
This makes efficient multi-mode QMs more realistic with 
current material properties, as we will discuss in the next 
Section. 



V. IMPLEMENTATION USING 
RARE-EARTH-ION-DOPED SOLIDS 

As storage material we here consider RE-doped crys- 
tals, which have inherently large inhomogeneous broad- 
enings (0.1 - 10 GHz) and very narrow homogenenous 
linewidths (0.1-100 kHz) at low temperatures (typically 
less than 4 K) 31]. Note that the typical homogeneous 
linewidths correspond to optical coherence times ranging 
from lis to ms. Still RF ensembles have relatively high 
absorption coefficients at normal doping concentrations 
(10-100 ppm), typically in the range 1-40 cm~^. RE crys- 
tals are therefore good candidates for the high resolution 
spectral shaping we are discussing here. Ref. [Slj con- 
tains an extensive review of optical properties of different 
RE-doped crystals and glasses. 

The creation of the periodic narrow structure is the 
crucial enabling step for the AFC QM. Narrow, absorbing 
peaks (< 50 kHz) have been isolated within the inhomo- 
geneous broadening by optical pumping (or spectral hole 
burning) [32, 33, 34, 35]. In aU stopped-light [22,,^] and 
CRIB experiments (S^, lU in RE solids to this date, a 
narrow absorption peak within a larger transmission hole 
has been created in order to have an approximate homo- 
geneous distribution of atoms. Briefly, a large transmis- 
sion hole is first created by spectral hole burning on the 
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optical transition while sweeping the laser frequency. The 
pumped atoms are stored in another hyperfine level with 
long population lifetime. A narrow ensemble of atoms 
is then created in the hole by repumping atoms from 
the storage state to_the initial state using a highly co- 
herent laser source [3d |. By repumping at different fre- 
quencies this method can be extended in order to cre- 
ate a series of peaks. One can also use other techniques 
in order to achieve frequency-selective optical pumping. 
For instance, one can use pairs of coherent pulses where 
each pair coherently transfers atoms to the excited state 
via a Ramsey interference [s^. Numerical simulations 
show that by accumulating many pairs, narrow, periodic 
structures can be created. This time-domain approach is 
interesting because in principle a time sequence can be 
created which has the appropriate frequency spectrum 
for creating an AFC with a certain finesse. Moreover 
a time-domain approach can be used for creating more 
complicated structures in frequency space, for instance 
multiple AFCs having different periodicity. 

The spin coherence times have the potential of be- 
ing significantly longer than the optical ones since it 
is free of decoherence due to spontaneous emission. In 
Pr'^+-doped ¥28105 crystals, one of the most studied 
RE-doped material, the spin coherence time is 500 /iS 
at zero magnetic field [38|, which can be extended to 
82 ms by applying an orientation-specific magnetic field 
[3^ . Fraval et al. [iO] has also shown that it is possible 
to use dynamic decoupling methods developed in NMR 
spectroscopy to further extend the coherence time to 30 
s. Yet another example is Eu'^+:Y2Si05 where 15 ms 
spin coherence time was measured at zero field [41], and 
Tm'^+iYAG where 300 /iS was obtained in a magnetic 
field Efficient transfer of population between the 

optically excited state and ground hyperfine states has 
been demonstrated in Pr'^+:Y2Si05, and manipulation 
of coherence between the optically excited state and the 
two spin ground states has recently been demonstrated 
in a state tomography experiment of an ensemble-based 
spin qubit in Pr'^+:Y2Si05 [43]. This shows the potential 
of spin states for long-term storage states for quantum 
information. 

We finally give a more detailed example based on 
Europium-doped ¥28105 crystals, to show the poten- 
tial for efficient multi-mode storage using an AFC based 
memory. Europium ions absorb at 580 nm and have an 
appropriate energy structure with three ground state lev- 
els and an optically excited state Q , as well as excellent 
optical and hyperfine ^41^ coherence times. Given 
the very narrow homogeneous linewidth in this material 
(122 Hz [4^), one could realistically create an AFC with 
7 = 2 kHz and A = 20 kHz (i.e. = 10). The to- 
tal AFC bandwidth we take to be 12 MHz, limited by 
hyperfine transition spacings, which results in Np — 600 
peaks. The maximum number of modes is then roughly 
N = Np/6 « 100. The absorption coefficient is about 
3-4 cm"-'^ [4^, so that 0? = 40 is within technical reach 
today using a multi-pass arrangement. Hence, one could 



then store 100 temporal modes in one single QM with an 
efficiency of 90%, which we verified by numerically simu- 
lating storage of 100 modes using the parameters above. 
We note that the efficiency stated above relates to the 
physics of the atomic frequency comb (see the model as- 
sumptions in sec. HI). A complete storage experiment 
including transfer to the spin state will be affected by 
the limitation in spin coherence time, for Eu'^+:Y2Si05 
j|spm_2g jjg^g been observed so far, and the trans- 

fer efficiency of the control fields (cf. discussion above). 

We conclude this section by observing that less studied 
RE-doped crystals could potentially store more modes 
due to larger hyperfine state separations, for instance 
Erbium-doped crystals absorbing at the telecommunica- 
tion wavelength 1.5 fim or Neodymium-doped crystals 
absorbing at 880 nm. 

VI. CONCLUSIONS 

We have presented a protocol for a quantum mem- 
ory that can store efficiently a large number of temporal 
modes in a single atomic ensemble by creating an atomic 
frequency comb. In materials with large optical inhomo- 
geneous broadenings and narrow homogeneous broaden- 
ings this allows for high multi-mode capacity at moderate 
optical absorption depths. We argue that RE-doped ma- 
terials have the necessary properties for implementing the 
proposed protocol. The inherent multi-mode property of 
the atomic frequency comb quantum memory promises 
to speed up the entanglement generation rates in long- 
distance quantum networks by several orders of magni- 
tude. 

The authors acknowledge useful discussions with Imam 
Usmani, Nicolas Sangouard, Thierry Chaneliere, Jean- 
Louis Le Gouet and Stefan KroU. This work was sup- 
ported by the Swiss NCCR Quantum Photonics and by 
the European Commission under the Integrated Project 
Qubit Applications (QAP). 

Note added. Since this paper was submitted for review, 
we have experimentally demonstrated a light-matter in- 
terface at the single photon level [i^ using an atomic 
frequency comb. 

APPENDIX A: ANALYTICAL CALCULATIONS 

1. Equations of motion 

In this appendix we present the analytical calcula- 
tions leading to the quantitative formulas in the Arti- 
cle. The basic equations describing the dynamics of the 
light-atoms system are the same as in section II of Ref. 
[29|. Both the field and the atomic polarization are de- 
scribed by (slowly varying) forward and backward modes, 
denoted Ef{z,t) and Eb{z,t) for the field and af{z,t;S) 
and ab{z,t;S) for the atomic polarization corresponding 
to atoms with resonance frequency 5. They satisfy the 
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following equations, cf. Eqs. (5a-b) and (6a-b) of Ref. 

faai. 



d_ 

dt 



(7f{z, t; S) = —iSaf{z, t; 5) + ipEf{z, t) (Al) 



ig-^+c-g-jEfiz,t)^zp I d5n{5)af{z,t-5) (A2) 



-Q^^b[z, t; S) = -iS(Tb{z, t; S) + ipEb{z, t) (A3) 

CO 

{-^^-c-^)Ebiz,t)^ip I d5n{5)ab{z,t-5). (A4) 



3. Absorption 

The absorption process is described by the forward 
equations (jAmA2p coupling Ef{z,t) and af{z,t;d). 
Plugging the solution of Eq. IjAip into Eq. (IA2|) one 
has the following equation for the field 



{Q-^+c—)Ef{z,t) = -pp j d6niS)e-^''x 

— CO 

t 

J dt'e'^'' Ef{z,t'). (A7) 



Here p is the dipole moment of the atomic transition, This can be rewritten as 
and p = glp, with go = \/^o/(2eol^), is the opti- 



cal frequency, and V is the quantization volume; n{5) is 

oo 

the atomic spectral distribution, with J dSn{S) = Na, 

where Na is the total number of atoms. See Ref. [291] for 
more details. Note that n{S) = NaG{6) and p = (i/Na 
in the notation of Ref. ^31 . 



2. Atomic spectral distribution and its Fourier 
transform 



Let us suppose that 



(A5) 



i.e. the atomic spectral distribution is a frequency comb 
with overall width F, peak separation A and individual 
peak width 7. Note that in the Article we talk about 
the Full- Width at Half-Maximum (FWHM) peak width 
7 = V8 In 27. 

We are interested in the regime where F ^ A 7, 
i.e. we have a significant number of well-separated peaks 
within the overall characteristic width. We define the 
finesse F of the atomic distribution as the ratio of the 
peak separation over the FW HM of the individual peak, 
which gives F = A/7 = A/(V81n27). 

In the following the Fourier transform fi{t) of n(5) will 
play an important role. One finds 



hit) oc e 



E 



(A6) 



which is also a series of peaks, with overall temporal 
width 1/7, peak separation 27r/A and individual peak 
width 1/F. Under the above assumptions, the peaks are 
again numerous (within the characteristic overall width) 
and well-separated. To summarize, both n{5) and fi(t) 
are pulse trains, with inversely proportional separations 
between the individual pulses, and the global shape of 
n{5) determines the local shape of fi{t) and vice versa. 



{gl+Cg-^)Ef{z,t)^-pp I dt'h{t-t')Ef{z,t') (A8) 



with 



hit) = / d6n{S)e 



-iSt 



(A9) 



Let us consider an incoming light pulse of temporal 
length T. Let us note right away that we are here consid- 
ering the regime r ^ L/c, where L is the length of the 
medium. This allows us to neglect temporal retardation 
effects, i.e. we can neglect the temporal derivative in Eqs. 
(jA7p and (jASp . The absorption process then has a du- 
ration that is essentially given by r. After this time the 
field is very close to zero, until the coherent photon-echo 
type re-emission at t = 27r/A. If we are interested only in 
the absorption process, we only need to consider values 
oi t ~ t' of order t in Eq. Assuming r < 27r/A, 

i.e. a pulse duration much shorter than the waiting time 
for the re-emission, only the central peak of h from Eq. 
[6l) will contribute, i.e. 



n(t-t')oce-(*-*')'rV2 



(AlO) 



for t—t' of order r. Assuming that the overall width of the 
spectral distribution F is much larger than the spectral 
width of the incoming pulse, i.e. Fr 3> 1, h{t — t') will 
act essentially like a delta function in the integral in Eq. 
[8|) . giving an approximate equation of the form 



— Ef{z,t) = --Ef{z,t), 



(All) 



with a constant absorption coefRcient a. This is easily 
solved, giving 



Efiz,t) ^ EfiO,t)e-'i^ 



(A12) 



Concerning the absolute strength of the absorption, one 
can show that under the above assumptions a = a ^^'^ ^ 
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where a is the absorption coefficient corresponding to the 
(central) peaks of n{5). This can be rewritten 



^ = 1.064^. 
4 In 2 F 



(A13) 



Hence, for the absorption process, the AFC is indis- 
tinguishable from a smooth distribution that would be 
obtained by a spectral averaging that does not resolve 
the individual peaks. Note that in the Article we fre- 
quently discuss the absorption depth d, here defined as 
d = aL. The effective absorption depth of the AFC 
is then d = aL, such that equation (|A13[) above corre- 
sponds to equation (4) in the Article. 



4. Re-emission 

The pair of tt pulses, which are used for storing the 
excitations in a long-lived ground state and then bring- 
ing them back to the excited state, also have the effect 
of changing the phase pattern of the atomic polarization 
from forward (c/) to backward (iTb), such that it now 
couples to the backward field Et, cf. Ref. [1^. We will 
see that this is essential for achieving potentially per- 
fect readout efficiency, in analogy with the case of CRIB 
memories [29j. We assume that the time spent in the 
ground state does not lead to any relative phases between 
the different frequency components (no inhomogeneous 
broadening in the ground state). From the point of view 
of the equations of motion, one can thus treat the two 
TT pulses as being applied essentially simultaneously. We 
will refer to the corresponding time as t = 0. Absorption 
thus occurred at negative times, and re-emission will oc- 
cur at positive times. The backward polarization created 
at t = is 



ab{z,0;S) = af(z,0;S) 



dt'e'^*'Ef{z,t'), (A14) 



where ab (ct/) is to be understood after (before) the two 
TT pulses. Eq. (|A14p gives the initial condition for Eqs. 
(3-4). We assume that the field Eb is zero at t = (the 
absorption process is completed, if there was a transmit- 
ted field it has left the medium) . This gives the following 
equation for the backward field 



{^^-c^)Eb{z,t)^-pp J dt'n{t-t')Ef{z,t') 

— oo 

t 

-pp J dt'n{t-t')Eb{z,t'). 



(A15) 

The first term on the right hand side is the source term 
which gives rise to the coherent re-emission, the second 



term describes the absorption of the backward field in 
analogy with the case treated in the previous section, 
see Eq. (|A8|) . The re-emission occurs because h{t — t') 
has a peak not only for t ^ t' around zero, but also for 
t — t' around 2tt/A. Note that the height of this peak 

is reduced by a factor e"^^-* ^ with respect to the 
central peak, cf. Eq. (|A6p . Around the corresponding 
point in time, one can apply arguments fully analogous 
to the ones in the previous subsection. Eq. (|A15p then 
simplifies to 



d_ 

dz 



Eb{z,t) 



'"Ef{z,t- 



2tt ^ a „ , 
^) + ^Ebiz,t). 



(A16) 

The Ef term has a factor a instead of ^ because the 
t' integral in Eq. (|A15p is over all of the effective delta 
function (given by the corresponding peak of n) for the 
Ef (source) term, whereas it is only over one half of the 
delta function for the Eb (absorption) term, as it was in 
the previous subsection for the absorption of Ef. 
The solution of Eq. (|A16P is 



Eb{z,t)^e-^"^"-^^"'^ j dz'aei^'''">Ef{z',t-^). 

(A17) 

Note that Et propagates backwards, so the output field 
is at z ~ 0. Using Eq. (|A12[) one has 

Eb{Q,t) ^ -Ef{0,t-^)e-^'(^->'/'x 



-EfiO,t^—)e 



(A18) 



It is interesting to note that the integral over z' has the 
same form as in the case of the (backward) CRIB memory 
[2^ . It corresponds to the coherent sum over all possible 
paths of the photon in the medium, where each path 
is characterized by its turning point z' . Even though 
there is absorption along the path, there is a constructive 
interference effect which leads to the output amplitude 
potentially being arbitrarily close in magnitude to the 
input amplitude (for large enough aL). 

This interference is not present if the output field is 
emitted in the forward direction, which would be the case 
if the two control fields were not applied. Then one would 
have an integral 

L 

J (i^'ae-""'/2e-"(^-"')/2 = aLe^"^^^, (A19) 



replacing the last factor in Eq. (jAlSp . If one then calcu- 
lates the efficiency, one finds that it can at most attain 
54% for 5L = 2, cf. [H. 
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For the backward propagating mode the achievable 
memory efficiency is given by, from Eq. (|A18[) . 



5. Phase due to frequency offset 



(A20) 



which can be rewritten in terms of the finesse F and peak 
absorption a of the atomic distribution n{S) as (cf. Eq. 
(5) in Article) 



rj = (1 - e V 4iT,2 )2g 
« (1 - e ~ye . 



(A21) 



There is a trade-off between absorption (favoring small 
F) and dephasing (favoring large F), however, for large 
enough aL and optimized F the efficiency can be made 
arbitrarily close to 100% (cf. Figure 2 and 3 in Article). 



Here we will comment on the phase of the output field 
as a function of the relative frequency detuning between 
the carrier frequency and the frequency of the peak in 
the AFC (any of them), cf. Eq. (5) in the Article. If the 
atomic distribution n(5) is replaced by n{5 — Aq), cor- 
responding to a small frequency offset, then the Fourier 
transform h{t) acquires a phase factor e*'^"*^ which leads 

to a supplementary phase factor e ~ for the re-emitted 
field, cf. Eq. (jA18|) . Note that this phase is modulus the 
peak separation A. 
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